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ff^ ' Abstract. In a recent paper | |astro-ph/0311364| ) Alam et al. argued that the SNe 

data of Tonry et al. 2003 and Barris et al. 2003 appear to favour DE which evolves in 
,,^^ ^ time, provided no other priors are invoked. (The effect of invoking priors such as the 

pg ' age of the Universe, the values of Hq and flom and CMB/LSS observations could modify 



% 



f~^ , this conclusion, as demonstrated in astro-ph/0403687j and other recent papers.) The 

>0 ' approach adopted by Alam et al. to reconstruct the properties of DE was severely (and, 

^S^ i as we shall show below - unfairly) criticized by Jonsson et al. in astro-ph/0404468 In 



■<!:;;j- I this paper we re-examine the parametrisation used in astro-ph/0311364 and show that, 

^D ■ contrary to the claims of Jonsson et al. , the results obtained from this reconstruction 

r^ , are robust and therefore representative of the true nature of dark energy. 

o 

^ ■ 1. Introduction 



The Universe appears to be accelerating and the nature of dark energy (DE), which 
drives this acceleration, is a subject of much current debate among cosmologists. The 
cosmological constant is the simplest possibility, but evolving dark energy models have 
also been suggested and one of the goals of current cosmological studies is to differentiate 
between these different types of models. Type la supernovae treated as standardized 
candles provided the first indications that the expansion of the universe is accelerating 
m 12] and one expects that the nature of dark energy will be further revealed by the 
study of these objects in conjunction with other data sets (CMB, LSS etc.). 

Recently there have been several data releases from the two supernova teams 
131111 El Ej, making the total number of supernovae almost double that known previously. 
Using these datasets, we reconstructed the dark energy density and equation of state 
in two papers [3 |H| (hereafter Papers I and II) in an attempt to understand the nature 
of dark energy. In these works, we found that the current supernova data appears 
to favour an evolving dark energy model with w < — 1 at present and that at 2a, 
the evolving model is at least as probable as the cosmological constant. This result 
soon found support in other works (see for instance [3 UHl UH 1121 UHl UH UH UHj etc.). 
Recently, in [T2j, doubts have been expressed as to the reliability of the fitting procedure 
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used in Papers I and II. Hence in this work we re-examine the reconstruction process 
used in these two papers to check whether the results obtained in Papers I & II were 
merely a consequence of the analysis or whether they do actually reflect a property of 
the universe. 



2. Reconstructing Dark Energy from Supernova Data 

The SNe observations measure the luminosity distance dilz). From this, information 
about the cosmological parameters may be obtained through the Hubble parameter, 
which in a spatially flat universe is related to the luminosity distance quite simply by 

[m 11111201 
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The dark energy density is defined as : 
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where poc = SHq/{8ttG) is the present day critical density of an FRW universe, and f^om 
is the present day matter density with respect to the critical density. 

Information extracted from SNe observations regarding diiz) therefore translates 
directly into knowledge of H{z), the dark energy density, and, through 
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into knowledge about the deceleration parameter of the universe and the equation of 
state of dark energy. 

One route to the reconstruction of DE lies in the construction of an ansatz for one 
of the three quantities: H{z),dL{z) or w{z). The ansatz must of course be sufficiently 
versatile to accommodate a large class of DE models. In Papers I and II we have used 
a polynomial fit to dark energy density of the form : 

H(x) 
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Note that this ansatz should not be considered as 



h{x) = 

where Aq + Ai + A2 

a truncated Taylor series for h'^{z). Rather, it is an interpolating fit for h'^{z) having 
the right behaviour for small and large values of z. The number of terms in this fit is 
sufficient given the amount and accuracy of the present supernovae data. With more 
and better data in the future, more terms with intermediate [e.g., , half- integer) powers 
of X may be added to it. 

This is equivalent to the following ansatz for DE density (with respect to the critical 
density) : 

Pbe{x) = Pde/poc = ^0 + Aix + A^x^ , (6) 
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Figure 1. Results from analysis of the 156 SNe "Gold" data from [0], using 
ansatz Q. ^Qm = 0.3 and h is marginalised over. Panel (a) shows the la, 2a, 2>a 
confidence levels in the (^1,^2) parameter space. The star marks the best-fit 
and the ACDM point is at the intersection of the dashed lines. The shaded 
grey region has /Ode < in the redshift range < 2; <, 2. Panel (b) shows 
the variation of dark energy density Pde(-z)/poc (where />oc = SHq/SttG is the 
present critical energy density) with redshift. Panel (c) shows the evolution of 
dark energy equation of state with redshift. In both the panels (b) and (c), the 
thick solid line represents the best-fit, the light grey contours represent the la 
confidence level, and the dark grey contours represent the 2a confidence levels. 
The dashed line in panels (b) and (c) represents ACDM, and the dotted line in 
panel (b) represents the matter density. 



which is exact for the cosmological constant w = —1 {Ai = A2 = 0) and for DE models 
with w = -2/3 {Ao = A2 = 0) and w = -1/3 {Aq = Ai = 0). 

The corresponding expression for the equation of state of DE is : 

Aix + 2A2X^ ,_, 

The likelihood for the parameters of the ansatz can be determined by minimising 
a x^-statistic: 



X^iHo,^Om,Pj) =J2 



[f^&t,i{Zi', Hq, fiom, Pj) — /^0,j 
^1 



(8) 



where /io,i = m^ — M = Slogrf^ + 25 is the extinction corrected distance modulus for 
SNe at redshift Zj, ctj is the uncertainty in the individual distance moduli (including the 
uncertainty in galaxy redshifts due to a peculiar velocity of 400 km/s), and pj are the 
parameters of the ansatz used. We assume a flat universe for our analysis but make no 
further assumptions on the nature of dark energy. We also marginalise over the nuisance 
parameter Hq by integrating the probability density e"^ ^^ over all values of Hq. 
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Figure 2. Results for analysis of SNe(Gold)+CMB data with ACDM-based 
WMAP priors of Oom = 0.27 ± 0.04 and h = 0.71 ± 0.06, using ansatz ©. 
Panel (a) shows the la, 2(t, 2>a confidence levels in the (yli, ^2) parameter space. 
The star marks the best-fit and the ACDM point is at the intersection of 
the dashed lines. The shaded grey region has pde < in the redshift range 
< z < 2. Panel (b) shows the variation of dark energy density Pde{z) / poc 
(where poc = SHq/SttG is the present critical energy density) with redshift. 
Panel (c) shows the evolution of dark energy equation of state with redshift. 
In both the panels (b) and (c), the thick solid line represents the best-fit, 
the light grey contours represent the la confidence level, and the dark grey 
contours represent the 2a confidence levels. The dashed line in panels (b) and 
(c) represents ACDM, and the dotted line in panel (b) represents the matter 
density. 



Figure [T] shows the results obtained by using the ansatz (0) for the "Gold" dataset 
of 156 SNe published in jB]. We find that the cosmological constant is just within the 
2(7 confidence level in the Ai — A2 plane. The variation of dark energy density and 
equation of state of dark energy is also shown, and both show indications of evolution. 
We note that these results are stable to the addition of extra terms such as x^ or x~^ 
terms to the ansatz. The region where the dark energy density becomes less than zero 
for the redshift range < z <^ 2 (the redshift range in which current supernova data is 
available) is shown as the shaded grey region in figure ^. We see that a portion of the 
3cr(99%) confidence level lies in the region where dark energy density goes to zero and 
therefore w blows up. 

Figure El shows the results of reconstruction when results from the WMAP 
experiment are incorporated and the priors Qq^ = 0.27 ± 0.04 and h = 0.71 ± 0.06 
are assumed. In this case the time-evolution of DE is considerably weaker and more in 
agreement with a cosmological constant. (Note however that evolving DE is preferred 
over a cosmological constant by the combined SNe+WMAP data if fiom — 0.4, h ^ 0.6 
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Figure 3. The fractional deviation Alog((ii^/7o)/log((iL//o) between actual 
value and that calculated using the ansatz (jSJ over redshift for different models 
of dark energy with i7om = 0.3. The solid lines represent quintessence tracker 
models for potential V = Vo/t/)", with a = 2 and 4. The dotted lines show 
the deviation for Chaplygin Gas models with k = 1 and 5 (where k is the 
ratio between CDM and Chaplygin gas densities at the commencement of 
the matter dominated epoch). The dot-dashed line represents the SUGRA 
potential, V = {M^+'^/^°') exp[i {(f)/Mpif], with M = 1.6 x 10-^Mpi,a = 11. 
The dashed horizontal line represents zero deviation from model values, which 
is true for ACDM and w = —l/3,w = —2/3 quiessence models. 

as shown in Paper II, 0.) 

3. Examining the reconstruction exercise for possible bias 

We now look at the reconstruction exercise in some detail to understand whether the 
results obtained in Papers I and II are biased due to inherent weaknesses in the ansatz 
used or whether these results can be trusted to reflect the true nature of dark energy. 



3.1. Reconstruction of different dark energy models 

In P7j , the ansatz (j3)) was criticised for being biased and giving misleading information 
about the nature of dark energy. We show here that this is not the case and the ansatz is 
sufliciently robust to be a reliable probe of the nature of dark energy. In the flgureEl we 
show the maximum deviation between the actual value of the luminosity distance and 
that calculated using the ansatz (jS)) for simulations of different physically well motivated 
models of dark energy such as quintessence, Chaplygin gas, and the SUGRA potential. 



Rejoinder to "No Evidence of Dark Energy Metamorphosis ", \astro-ph/04 044 68 6 



From (0) it is clear that for w = —1, w = —2/3, and w = —1/3, this ansatz returns 
exact values. However, even for models for which it is not exact, the ansatz (jSj) still 
recovers the luminosity distance to within 0.5% accuracy in the redshift range relevant 
for SNe observations. These results reassure us that for a large class of models, the 
ansatz can be trusted to recover cosmological quantities to a high degree of accuracy. 
(Other tests of the ansatz (J3J have been reported in J221I231IZ1-) 

3.2. Diverging equation of state 

In P2! the ansatz © has been criticized on the grounds that it allows the DE density to 
become negative during the course of cosmological evolution. We would like to emphasise 
that this is not necessarily a bad thing since there is no a-priori reason for the DE density 
to remain positive throughout its evolution and models do exist in which p^E < for a 
substantial duration of time - see for instance [2111231231221121] etc. Indeed the purpose 
of the construction of an ansatz such as ^ is to endow it with sufficient flexibility so 
that it is able to reproduce the behaviour of a sufficiently large class of DE models. | In 
addition, a negative value for the dark energy density could also be an indication that 
the matter density has been over estimated, so one should not arbitrarily exclude the 
region of parameter space where pde < as suggested by [l2|. Hence we have not a 
priori put a constraint of positivity on the dark energy density in our analysis. However, 
one point to note is that the region where dark energy density vanishes is very close to 
the cosmological constant {Ai = A2 = 0) in parameter space, so there is a chance that 
if the best-fit to the data is near this point then the best-fit w or errors in it may blow 
up. This need not always happen, for instance in the figure |21 the best-fit is close to 
the cosmological constant, but the errors on w are well-behaved. The same is seen in 
figure 6 of Paper I. However, in case the best-fit equation of state or the errors in it do 
blow up, we naturally have to be careful about whether this truly indicates some exotic 
form of dark energy, or whether it is more consistent with the cosmological constant. It 
is also important to point out that the equation of state parameter need not always be 
the best characteristic of dark energy |^ , and in such cases it is useful to characterize 
DE using other cosmological parameters (which do not exhibit this divergence), such 
as the deceleration parameter and the Statefinder pair |221 1231 • In this context it is 
worth noting that the pressure of dark energy : P/poc = — Pde + Aix/3 + 2^2X^/3 is 
well-behaved for this ansatz even when pde < 0, and this quantity can therefore be used 
to draw conclusions on the nature of dark energy. 

3.3. Non-monotonic errors on the equation of state 

The errors on the equation of state in figure^c) appear to be non-monotonic. This seems 
somewhat counter-intuitive. Low redshift behaviour of the equation of state affects the 

I It is clear that no finite parameter ansatz will will ever be able to reproduce all possible behaviour of 
DE. For instance the ansatz (O is unlikely to give good results if applied to DE which evolves extremely 
rapidly with time such as the rapidly oscillating models of |29| . 
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Figure 4. The variation of the equation of state of dark energy w{z) over 
redshift for the ansatz ©. The thick sohd hne shows the best-fit and the hght 
grey contour represents the lo" confidence level around the best-fit. The dashed 
horizontal line denotes ACDM. fiom = 0.3 is assumed. Errors are calculated by 
the Fisher matrix approach, using Eg 1121 



luminosity distance at all higher redshifts, while high redshift behaviour affects fewer 
points. This leads to an expectation that high- 2; behaviour of the equation of state 
should be poorly constrained as opposed to the low-2 behaviour. However, we saw in 
Papers I & II that the errors in w{z) actually decrease with redshift. We investigate 
this matter using the Fisher matrix error bars. 

In an analysis which uses an ansatz with n parameters pj, the Fisher information 
matrix is defined to be 

^'' ^ \d^) ' ^^^ 

where L = —logC, C being the likelihood . For an unbiased estimator, the errors on the 
parameters will follow the Cramer- Rao inequality : Api > I/^/FH. 

Since the likelihood function is approximately Gaussian near the maximum 
likelihood (ML) point, the covariance matrix for a maximum likelihood estimator is 
given by 

)ij 
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The Fisher information matrix is therefore simply the expectation value of the inverse 
of the covariance matrix at the ML-point. 

Given the covariance matrix, the error on any cosmological quantity Q{pi) is given 
by: 
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Thus the nature of the errors on a quantity will depend essentially on the manner in 
which it is related to the parameters of the system. 
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The errors on the equation of state for the ansatz (0) can now be calculated using 
equations (0) and (fTTj) . and has the somewhat complicated expression : 



2 . , _ X^iflC, + 2/1/2C12 + /IC22] 



-^ ^ 9[1 - fiom + Ai(x - 1) + A2(X2 - 1)]4 ' ^ ^ 



where 



/l = l-fi0m-A2(x-l)^ 

/2 = 2x(l-fiom)+^l(x-l)^ 

and Ai, A2 are the best-fit values of the parameters. 

From this expression it is difficult to predict whether the errors should increase 
with redshift or not. Indeed that would depend on the value and sign of the quantities 
Cii. We calculate the errors around the best-fit for the 156 SNe "Gold" dataset of 
using this procedure that the la confidence levels thus calculated on the equation of 
state (shown in figure |^ are almost identical to those shown in figure Q^c). § 

Thus the ansatz (0) need not show monotonically increasing errors in w with redshift 
and in this way may differ from other fitting functions, for instance approximations to 
the equation of state of dark energy as proposed in ^^ or ^^, which may show just 
the opposite tendency. This indicates that the nature of error bars are affected by 
the quantity being approximated. In the limit of infinite terms in the expansion of 
various quantities all the methods should produce identical result. The practical need 
for truncating these expansions make these approximations slightly different from each 
other. More specifically, we require setting of priors 



f{z)=Y.(^^^^ (13) 

i=0 

an =0; {n>Np) (14) 

where f{z) could be H{z), w{z) or any other physical quantity and Np is the chosen 
number of parameters. The non-linear priors in the above equation make different finite 
expansions inequivalent. Since we do not know for certain if the underlying model for the 
accelerating expansion involves an energy component with negative pressure in a FRW 
setting we are forced to choose one of the alternatives for approximations. We hope 
that with increasingly high quality data the effect of such truncations will eventually 
disappear. 

3.4- Stability of the dark energy equation of state 

The dark energy equation of state for the ansatz ^ can have many different forms 
depending on the values of Ai,A2. This ansatz gives us an exact representation of the 
cosmological constant {Ai = A2 = 0) and it is an unbiased estimator for the ACDM 
model. However, it is also true that the statistical noise inherently present in the data 

§ One should note that the Fisher matrix approach works best for the cases where the hkelihood is 
symmetric aromid the best-fit. For asymmetric hkcHhood curves, this method should not be used. 
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will never allow us to measure Ai, A2 to be exactly equal to zero even in this case. Thus 
if the underlying model of the universe is the cosmo logical constant, then by fitting an 
ansatz which has more parameters than required, we are fitting a bit of noise as well. 
This fact is true for any ansatz and not just (0). The deviation of the best-fit from the 
cosmological constant in such a case is given by : 

i.e. the w\ curve would increasingly depart from the cosmological constant with redshift. 
The same tendency would be seen if we consider a linear ansatz in the equation of state : 
w = wq+wiz. This would also diverge from the cosmological constant {wq = —1, wi = 0) 
at high redshifts for even a small non-zero 6wi, since 6wa ~ 6wo + 6wiZ. Such 
divergence is the inevitable price we pay for not knowing the "true" nature of dark 
energy and for using an ansatz to explore it. In practice one is free to choose either 
of two possibilities: (i) assume that the underlying model has certain characteristics 
(such as constant equation of state) and thus face the danger of losing information- the 
pitfalls of using this approach were highlighted in P2|, or (ii) we may choose to make no 
such assumptions and consider models with a somewhat larger number of parameters. 
We have adopted (ii) in Papers I & II and believe that, given the redshift range and 
statistical noise of the data currently available, it is not unreasonable to assume that 
the errors for the ansatz © will not be large enough to cause the best-fit to diverge 
significantly from the underlying model irrespective of whether the model is that of the 
cosmological constant or some other form of dark energy. This is demonstrated in figure 
[7|^a) for simulations of the cosmological constant model. 

3.5. Correlation of the parameters Ai, A2 

In our analysis of the supernova data we find that for a fixed value of fiom, the parameters 
Ai, A2 which are favoured by the data are strongly correlated. For Qq^ = 0.3, from the 
confidence levels of figure ^a) we may conclude that approximately, Ai ~ —2.AA2 for 
points favoured by the data. Contrary to what is claimed in ^\ this correlation is not 
something that has been put into the analysis by hand, rather it is a consequence of the 
observational data. Other fits, such as the hinder fit to the equation of state [SO] , have 
also shown a strong correlation between its parameters. We need to examine what effect, 
if any, this correlation has on the nature of dark energy. From expression ((Zj) for the 
equation of state for dark energy, we see that w = — 1 at a redshift oi z = —A1/2A2 — 1. 
Therefore, for r^om = 0.3, w = —1 at 2; ~ 0.2. This implies that since most curves 
would cross w = —1 at approximately this redshift, there will be a "sweet-spot" at 
z ~ 0.2 for the equation of state. This can be clearly seen in the figure ^c). The 
presence of sweet-spots at different redshifts for different ansatz is something which has 
been studied earlier, for example in [32], and the community agrees that the presence or 
absence of a sweet-spot is simply a consequence of the ansatz used and does not signify 
any extra information about the underlying model. Does the fact that w crosses —1 at 
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Figure 5. The variation of the equation of state with redshift for different 
values of the parameters Ai,A2 satisfying Ai = —2.4^2. The thick dashed hue 
denotes ACDM. 



this redshift mean that the ansatz forces the equation of state to evolve in a particular 
direction? In figure El we plot the behaviour of w for different values of the parameters 
Ai, A2 which follow the constraint Ai = —2AA2. We see that the cosmological constant 
is allowed by this constraint, as are many other models. One can have extremely mild 
evolution of w, or more rapid evolution of w, one can have w which blows up, w which 
goes from a more negative (< — 1) value to a less negative (> — 1) value or vice versa. 
In short this constraint does not force the equation of state to evolve at a particular 
rate in a particular direction, but leaves a great deal of freedom. Also, it is important 
to note that this is not a constraint that has been imposed by hand in the analysis, it 
is simply an outcome of the analysis. There are many points in the Ai,A2 parameter 
space for which the equation of state is always > — 1 and shows little evolution. The 
ansatz does not preclude these points, it is the data which chooses a particular evolving 
dark energy model over other possibilities. We therefore conclude that the correlation 
seen in the parameters Ai, A2 does not signify any bias in the results and (contrary to 
the claims of J7]) the ansatz does not "force" evolution of the dark energy equation of 
state. 



3. 6. Confirmation from other fitting functions 

In Paper I, we showed that the results obtained by the ansatz (0) can also be reproduced 
using other fitting functions, such as the fits to the equation of state used in [20] and 
PT] . For both these fits, the equation of state shows an evolution which is remarkably 
similar to that found using ©. It is interesting that our results have found support in 
the independent analysis carried out by other groups which also show that evolving DE 
provides a fit that is as good (or better) than that provided by a cosmological constant 
to the current supernova data [01 CHI HH 1121 CHI- Therefore there seems to be a consensus 
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among the different groups that, based on the SNe data alone, an evolving dark energy 
model is a plausible alternative to the cosmo logical constant. (As emphasised earlier in 
this paper and in Papers I and II, the properties of dark energy depend crucially upon the 
manner in which the SNe data is sampled and also upon which other sets of observations 
are used in conjunction with type la supernovae. For instance, the evolution of DE is 
least if one uses SNe+WMAP results together with the priors VLq^ = 0.3, h = 0.7; see 
also ^.) 

4. Simulation of a ACDM model 

To test the reliability of our ansatz, we attempt to reconstruct a fiducial cosmology 
from simulated data. We assume that the background model is a fiat ACDM universe 
with fiom = 0.3, and generate 1000 datasets, each consisting of 156 SNe with the same 
redshift distribution and magnitude errors as those in the "Gold" dataset of Ej. We 
first analyse these datasets using the polynomial fit to dark energy density, Eq (jSJ. The 
best-fit values of Ai,A2 for these 1000 sets are shown in figure IHl The mean of this 
set of 1000 points is very close to the cosmological constant: Ai = —0.01,^2 = 0.005 
{Ai = A2 = for ACDM). The best-fit points consistent with the cosmological constant 
at 68% confidence level are indicated by filled circles, points consistent at 95% confidence 
level are shown by crosses, those consistent at 99% confidence level are shown by filled 
triangles, and those inconsistent at 99% confidence level are shown by the open circles. 
By determining confidence levels around the best-fit for each individual dataset we find 
that about 75% of the best-fit points are consistent with the cosmological constant at 
Icr. We also note that some of the best-fit points are in the region where pde < 0. We 
now choose three of these 1000 datasets and calculate pde and w for the corresponding 
cosmology. The three best-fit points chosen are shown- the first (hereafter called dataset 
(a)) is consistent within 68% confidence level, the next (dataset (b)) is consistent at 99% 
confidence level and the last (dataset (c)) is inconsistent at 99% confidence level with 
the cosmological constant. The results are shown in figures [7| and |H1 We find that for the 
dataset (a), pde and w are consistent with the cosmological constant, even though the 
errors in w blow up. For the datasets (b) and (c), the nature of dark energy density and 
dark energy equation of state appear to be very different from that of the cosmological 
constant. Such an example was used in figure 4 of f^j to make the claim that the 
ansatz (0) leads us to conclude that w is evolving even though the fiducial model is that 
of the cosmological constant. 

Let us now see if this is really correct and whether our results were indeed a 
construct of the ansatz used. 

• Let us consider the same fiducial model (ACDM) but an alternative ansatz to 
reconstruct dark energy. As our first example we consider a cosmological constant 
model but make no assumption about the flatness of the universe. The luminosity 
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Figure 6. A scatter plot of the best-fit values Ai,A2 of the polynomial fit 
to DE density, Eq (jSJ, for 1000 simulated data sets in a flat universe with 
a cosmological constant. Oom = 0.3 is assumed. The filled circles represent 
simulated datasets consistent with ACDM at the 68% confidence level, the 
crosses represents those consistent at the 95% confidence level, the filled 
triangles are consistent at the 99% confidence level, while the open circles are 
inconsistent at the 99% confidence level. Three of these datasets used in further 
analysis are shown by (a) a grey circle, (b) a grey triangle, and (c) a grey open 
circle. The grey shaded region has pde < in the redshift range < z < 2. 
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Figure 7. The variation of dark energy density /9de(-2)//00c (where poc = 
SHq/SttG is the present critical energy density) with redshift, using three sets 
of simulated data for a fiducial flat ACDM model. The polynomial fit to DE 
density, Eq © with fiom = 0.3 in a flat universe has been used. In each panel, 
the thick solid line shows the best-fit, the light grey contour represents the la 
confidence level, and the dark grey contour represents the 2a confidence level 
around the best-fit. The dotted line denotes matter density r2om(l + -z)^, and 
the dashed horizontal line denotes ACDM. 
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Figure 8. The variation of equation of state of dark energy w{z) with redshift 
for r^om = 0.3 in a fiat universe using three sets of simulated data for the 
polynomial fit to DE density, Eq ©. The background model is fiat ACDM. 
In each panel, the thick solid line shows the best-fit, the light grey contour 
represents the la confidence level, and the dark grey contour represents the 
2(T confidence level around the best-fit. The dashed horizontal line denotes 
ACDM. 
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Figure 9. The la, 2a, 3a confidence levels in fionn^A for the three sets of 
simulated data (with a flat cosmological constant background model, J7om + 
Oa = 1) using a cosmological constant model, Eq (|16|) for the analysis. In 
each panel, a star marks the best-fit. The intersection of the dashed lines 
denotes ACDM. The solid line denotes the flat universe with Oom + ^A = 1- 
All universes above this line are closed and those below it are open. 



distance in this case is given by : 

, , , c{l+z) 
dL[z) 




dz 



(16) 



Ho^^W [' 'JO ^{i + zy{i + no^z)-zi2 + z)nA 

(Qf^ = 1 — fiom — ^a; <S = sin, 1, sinh for Qf, <, =, > 0) . 

The parameters of the system are fiom, ^a- The results are shown in the figure IHl 
We see that in case of dataset (a), the errors on Oom, ^a are completely consistent 
with the cosmological constant in a fiat universe, but for both the datasets (b) and 
(c) a closed model is the preferred model and the fiat model (which is the true 
background model) is excluded at la and 2cr respectively ! 
We may also consider different fitting functions for w in a fiat universe. 
Consider the popular linear fitting function for the equation of state : 



w{z) = Wq + WiZ (17) 

H\Z) = ifo'[^Om(l + Zf + (1 - fiom)(l + ;^)3(l+-0-l)e3^i.] . 

This is the ansatz used in |2j to parametrise dark energy. The resultant confidence 
levels in Wo,Wi are shown in the figure HUl Once again, we see that only in case 
of the dataset (a) does the parametrisation come close to reconstructing the actual 
model. For the dataset (b), the cosmological constant is ruled out at la and for 
the dataset (c) it is ruled out at 2a. 
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Figure 10. The la, 2a, 3a confidence levels in wq,wi using three sets of 
simulated data for the linear fit to equation of state of DE, Eq ()17() . r^om = 0.3 
and a flat universe are assumed. Background model is flat ACDM. In each 
panel, a star marks the best-fit. The intersection of the dashed lines denotes 
ACDM. 
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Figure 11. The variation of equation of state of dark energy w{z) with redshift 
for J7om = 0.3 and a flat universe using three sets of simulated data for Linder's 
fit to the equation of state of DE, Eq ()18|) . Background model is flat ACDM. 
In each panel, the thick solid line shows the best-fit, the light grey contour 
represents the la confidence level, and the dark grey contour represents the 
2a confidence level around the best-fit. The dashed horizontal line denotes 
ACDM. 
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• Finally consider the parametrisation suggested in [HO] : 

w{z) =wo + ^ (18) 

H\z) = ifo'[^Om(l + Zf + (1 - fiom)(l + z)3(l+-0+-l)e-3«.i./(l+.)] _ 

The evolution of the equation of state with redshift in this case for the three datasets 
is shown in the figure ITTl We see that in the first panel (a), the evolution is very 
moderate and is entirely consistent with a cosmological constant. In the second and 
third panels however, we see significant evolution of the equation of state, which is 
consistent with what is seen in the panels (b) and (c) of figure HI and inconsistent 
with a cosmological constant model. 

The above examples clearly indicate that the different fits give results which are 
mutually consistent and that, contrary to the claims made in ^2], the ansatz (jSJ does not 
"force" evolution of the equation of state of dark energy when the data is commensurate 
with the cosmological constant. If the datasets for which ^ shows evolution of DE are 
analysed using other fits, similar (misleading) conclusions about the behaviour of DE 
are reached. So these results are dependent not so much on the ansatz as on the 
dataset itself. Therefore it is in principle possible that even with a fiat ACDM fiducial 
model having errors commensurate with current observations one can reconstruct dark 
energy which evolves, or a universe that is closed ! However, one must remember that 
out of the 1000 realisations of such a model, about 75% are broadly consistent with 
the cosmological constant model for the ansatz ^ and if these datasets are analysed 
then the correct model will be recovered within la, as in case of the dataset (a). The 
datasets (b) and (c) which give rise to incorrect results are much less likely to occur. As 
the quality of data improves, the likelihood of obtaining datasets like (b) and (c) will 
further decrease. Given the current observational standards, it is reasonable to assume 
that the actual data which we have is consistent with the underlying model and results 
obtained from it are reliable. As more data becomes available, we will be able to discern 
the nature of the universe with greater accuracy. 

5. Conclusion 

We have shown that the results of Papers I and II are robust and that the conclusions 
drawn there are not due to any bias in the ansatz used. Contrary to the claims in ^\ 
the ansatz (jSJ does not inherently "force" or "favour" evolving dark energy models over 
the cosmological constant. This ansatz can reproduce the behaviour of a large class of 
DE models fairly accurately and as such it can be used with confidence to predict the 
nature of the universe using currently available data. The conclusion that the current 
supernova data favours the evolving dark energy models over the cosmological constant 
at la still holds and this result has found support in the work of other independent 
groups mH HH 1121 Uni HH Uni Uni when no other data sets are used. Better quality 
data expected in the future from different cosmology experiments (SNe, CMB, LSS etc.) 
will allow us to draw firmer conclusions about the nature of dark energy. 
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